ALGORITHM IN LATTICE C TO GENERATE S(n) 


by Jim Duncan 
9 Ryeground Lane 
Formby 
Liverpool 
L37 TEG ENGLAND 


Computer: Atari 1040ST 
Run time: ca. 50 hrs to generate S(n) for 1 000 000 numbers 


#include <stdio.h> 
Hinclude <math.h*? 


unsigned long int pst,ast,s,t; 


main() 
{ 


} 


long int n; 
FILE *fp; 


printf ("input n (1 <= n < 2 147 483 647)\n"); 
scanf("Zld", &n);3 

fp = fopen("PRN:“, "w"); 

fprintf(fp, "n = Zid Sín) = ZAld\n", n, smaran(n))s 
fclose(fp}; 


smar an (m) 
unsigned long int am; 


€ 


unsigned long int mst,p,a,fact; 
double r; 


if (m == 1) 
return(0); 
else ¢ /* STANDARD FORM of m #/ 
r= mst = m; 
p = 1; 
fact = 1; 
while ( ++p <= sqrt(r)) í 
a = Q0; 
while (mst 4 p == 0) ¢ 
mst = mst/p; 
att; 
+ 
r = mst; 
if (a > 0) í< 
pst = ps 
ast = a; 
t = s = O; 


11 


3 


tors! () 
t 


12 


tors2() 


tors3() 


{ 


/* 
/%* 


tors3(); /* find smallest factorial (t) 


/* p^a divsor #/ 
if (t > fact) 
fact = t; 
2 

2 
if (mst > fact) 

fact = mst; 
return(fact); 
2 


unsigned long int i; 


test number is pst^ast #/ 
s is the difference between a factorial number t and ast#{pst-1) #/ 


with #/ 


/* s forms a pattern which determines the smallest value of t for which */ 


/* 


the test number is a divisor #/ 


i = Q; 
while (++i < pstæpst && t-s < ast#¥(pst-1)) í 
if (1 X% pst == 0) 
s = s-pst+2; 
else 
s+t+; 
t += pst; 
2 


unsigned long int i; 


torsi(); 
i = 0; 
while (++i < pstæpst && t-s < ast#(pst-1)) í< 
if (i % pst == 0) 
s = s—S*pst+43 
else 
S$ = s~2epst+3; 
t += pst; 
torsi()s 
} 


unsigned long int i; 


tors2(); 
i = 0; 
while (++i < pstæpst && t-s < ast#(pst-1)) í 
if (i X pst == 0) 
s = s-S*pst+d6; 
else 
s = s-4#pstt+5;5 
t += pst; 
tors2(); 
2 


